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Abstract 

The vacuum energy density (Casimir energy) corresponding to a massless scalar quantum 
field living in different universes (mainly no-boundary ones), in several dimensions, is cal- 
culated. Hawking's zeta function regularization procedure supplemented with a very simple 
binomial expansion is shown to be a rigorous and well suited method for performing the 
analysis. It is compared with other, much more involved techniques. The principal-part 
prescription is used to deal with the poles that eventually appear. Results of the analysis 
are the absence of poles at four dimensions (for a 4d Riemann sphere and for a 4d cylinder of 
3d Riemann spherical section), the total coincidence of the results corresponding to a 3d and 
a 4d cylinder (the first after pole subtraction), and the fact that the vacuum energy density 
for cylinders is (in absolute value) over an order of magnitude smaller than for spheres of 
the same dimension. 
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1 Introduction 



The investigation of the vacuum energy density (or Casimir energy) |J corresponding to a 
quantum field theory defined in a certain manifold (spacetime) with or without boundary, 
is one of the most basic issues of quantum field theory. The Casimir force is probably the 
simplest and most spectacular of the different manifestations of the vacuum energy — better, 
of the modification of the vacuum energy when boundaries, curvature or a background field 
are superimposed to the manifold jl[ To understand the correspondence between the 
sign of the effect (associated with an attractive or repulsive force, depending on the sign 
being negative or positive, respectively) and the specific topology of the boundary is a most 
challenging point which requires explanation. A historical failure in relation with this issue 
ocurred some years ago, when the Casimir force was thought to contain the explanation of 
quark confinement in the context of the bag model. The idea was that this vacuum effect 
would provide the attractive force giving tension to the bag that contained the quarks. But it 
turns out that the calculation of the force for curved boundaries is very tricky, plagued with 
infinities and needing apropriate regularization. Sometimes cut-offs remain, zeta-function 
regularization is said to be insufficient, and it is difficult to extract uncontroverted, physically 
meaningful results ||. Actually, this subject has been the source of several sound errors in 
the scientific literature. In any case, when it was proven, without doubt, that the force for 
a closed sphere in three-dimensional space is repulsive (as it is for a closed cube, and not 
attractive as for a pair of plates) the brilliant idea of supplementing the bag model with the 
Casimir force had to be abandoned. 

Notwistanding that, the presence of the Casimir force in very different phenomena of 
condense matter, solid state and laser physics has been rigorously stablished, both theoreti- 
cally and experimentally ||. Its relevance for possible models of our universe is also without 
discussion. Maybe owing to the difficulties encountered when trying to give a plausible an- 
swer to the question: what are the boundary conditions of our universe ?, some of the most 
popular models of spacetime nowadays are given by manifolds without boundaries. 

Riemann spheres are to be counted among the simplest and most important of these 
manifolds. This is the reason why we have chosen them as the spaces of our models. We will 
study these models in different dimensions, with the hope to find some characteristic that 
may singularize some of the manifolds considered and some particular dimension, among all 
of them. 

From a more technical point of view, the specifications of our study are as follows. We 
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shall use the most elegant, rigorous and simple of the regularizations known to date, namely 
zeta function regularization |3]]. Moreover, we are going to supplement it with a most easy 
technique, which is binomial expansion [BJ. Then we shall compare our findings with those 
coming from other approaches, which turn out to be much more artificial, lengthy and less 
rigorous. Hard to believe as this statement may seem, it is the plain truth. Little merit about 
it, again this is just a confirmation of the famous general principle, atributed to Einstein, 
which asserts that nature always follows the most simple path among the ones that are 
available to her. That such path turns out to be in this case, at the same time, the most 
rigorous mathematically is however a rewarding small surprise (since sometimes physical 
intuition and beauty has been associated with mathematical sloppiness). Even from the 
pure calculational point of view, the method here developed yields very rapidly convergent 
expressions which allow us to obtain 6 or 8 digit precission with just a few first terms of 
a series, in a home computer and with a standard computation package (Mathematica, for 
instance, but never use it for dealing with the Hurwitz zeta function ((s,a) when a > 1). 
As side products of our analysis, some asymptotic expansions for Hurwitz and Epstein zeta 
functions, that had been obtained by the author previously, are here chequed and numerically 
contrasted with other results in the literature. As has been pointed out already, this is a 
necessary exercise in this field, because of the discrepancies and errors that so frequently 
appear. In particular, a small table of derivatives of the Riemann zeta function, which are 
repeatedly used in our numerical calculations, is given. 

In Sect. 2 we summarize our method of zeta function regularization and describe the 
way in which the binomial expansion is used. We compare it with other approaches that 
have been employed, in other to prove the advantages of our procedure. In Sect. 3 we 
obtain the Casimir energy density for Riemann spheres in d = 1,2,3,4 dimensions. These 
are manifolds without boundary, but also the situation when one has a half such manifold 
with Dirichlet and Neumann boundary conditions, respectively, is considered for comparison. 
Cylindrical manifolds whose sections are Riemann spheres, with total dimension d = 2, 3, 4, 5 
are investigated in Sect. 4. The discussion of our results and the conclusions of our analysis 
are given in Sect. 5. 
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2 Analytical approach to the zeta functions for Rie- 
mann surfaces 



As has been observed elsewhere ||, the formula J?) (see also ||) 
F(s;a,b) ee J2[(n + a) 2 + b] = ^E l j , J ^"C(-2n, a) (1) 

n=l 1 V S / n=0 n - 

v^r( S - l/2) 1/2 _, 2vr* 1/4 _ 5/2 ~ n i/2-s K j 2 ^ 

2T(s) r(s) 1 ; ^ 7 1 ; 

and its generalizations to higher dimensional manifolds J7|, are extremely useful (owing to 
exponentially quick convergence) when b > 0. However, it turns out to be rather difficult to 
apply them when b < 0. And this is usually the relevant case which appears when one tries 
to calculate, e.g., the determinants of the Laplacian operators on Riemann spheres by using 
zeta functions. One has to regularize (i.e., analytically continue) expressions of the general 
form 

oo 

f(s;a,b,c) = Y, l ~ S+b ( l + *r S+C ' (2) 
i=i 

where a turns out to be positive. 

An alternative expression to (Q) has been obtained in ||, by making use, essentially, 
of the usual integral representation for the Hurwitz zeta function together with the Mellin 
transformation (actually the same ingredients, aside from adequate series commutation, 
which were used in [|7j for the derivation of @). Specifically, for 

oo 

C (n) (s) = E [l~ s (l + 2n + I) 1 ' 8 + ^{l + 2n + l)~ s ] , (3) 
l=i 

the following interesting result was obtained 



as 



i 2n+l 

= 4C'(-1) - ~(2n + l) 2 + E (2k - 2n - 1) log*, (4) 

s=0 Z k=l 



from which, in particular, for the zeta function of the Laplacian on the hemisphere with 
Dirichlet and Neumann boundary conditions, respectively, 



one gets 



Cd(s) = e m + i)r> us) = eg + my + i)r, (5) 
i=i 1=1 



C D (0) = 2C'(-1) + I log(2vr) - I C N (0) = 2C(-1) - 5 log(27r) - J. (6) 
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Those are nice results, indeed. However, the derivation of these expressions in || is not 
free from difficulties. In fact, to start with, the analysis is rather lengthy and, on the other 
hand, a highly arbitrary, additional regularization is needed at some point: a certain infrared 
convergence factor t s (t is the integration variable) must be introduced and the exponent s 
can be let to go to zero only after performing a convenient combination of different terms. 
That these manipulations are not so obviously accepted (even if in the end they turn out to be 
right, as is the case here) is proven by the continuous recourse to specific checking of the final 
numbers with well-known results ||, 110, O. In certain cases, in fact, use of these integral 
transforms can lead to discrepancies with already known results. So, the value of ('(—1) 
obtained in flTJ] as a byproduct of an original method there developed seems to be quite 
far from the best accepted value as given, for instance, in ]13| (see expressions fl50f ) below). 
Our value is coincident (at least to 8 digits) with the one obtained by Salomonson using a 
different formula ||14|| . On the other hand, however, the remarkable expression obtained in 
12| for £'(— 2) appears to be right (see ([j0|)). 

A mathematically clean, rigorous and, at the same time, much more simple procedure 
to deal with any expression of the general form (||) can be devised which makes use of the 



simplest of ideas: binomial expansion ||. This goes as follows 

oo 

f(s;a,b,c) = ]Tr^+ c (/ + ar+ c = ; 

i=i 

OO OO TV1 r> I ^ 

-a k l~ 2s ~ k+h+c 



i=i 



- 2s+b {i + a r l )- s+c 

1-s + c) 



£{}+£{} 

1=1 l=[a]+l 



(7) 



being [a] the integer part of a, so that g(s) is an integer function of s, while the second, 
truncated series is absolutely convergent (since al' 1 < 1 there). The final result is 

m 



X 



C(2s + k-b-c) -J^l 



-2s-k+b+c 



1=1 



(8) 



In particular, 



m + c) 



/(o; a ,6,c) = E''(i^r'r + g Hr(1 _ fc + c) 



l=i 



C(k-b-c)-j2i 



-k+b+c 



(9) 



and 



f'(0-a,b,c) = -f^togZ + logCl + oT 1 )]^! 



i=i 



at 



-l\c 



+ E 

fc=0 



V>(1 - k + c) --0(6) 



+ 



r(i - k + c 

2T(l + c) 



C(k-b-c)-J2l 



1=1 



ra-k + c) 



('(k-b-c) + ]Tr fc+b+c logZ 



(10) 



In fact this last formula is a bit tricky, and has to be modified (in general) when b and c are 
integers. Then, the derivative for the particular value k = b + 1 must be done with special 
care, by performing first expansions around the poles and zeros of these functions at s = 
(this will be shown in full detail below). 

For the zeta function (d(s), eq. (H), we obtain 

T(l-s) 



<n(s) = j:i 1 - 2s (i + r 1 ) 

1=1 



2- s + y 

k\T(l - s - k) 



[C(2s + &-!)-!], (11) 



and 



&(0)=2C'(-l) + | + |-lo 



g2 -Y(-i) k C(A;) -. 



(12) 



The power of our simple method will be now demonstrated by showing that this last expres- 
sion coincides with the first of (|5]). In fact, we have 



E(-i) fe 



k=2 



k+l ^ 2 (k + l)r(fc) Jo 



oo + 

dt 



k-l 



e* - 1 



with 



& 



k=2 



It is easy to see that ip'(u) = u[ip(l) — ip(u + 1)] and integrating, ip(l) 
(1/2) log(27r). Finally, adding the rest of the series: 



E 

k=2 



(-1) 1 , 

ITT--2 +log2 ' 



(13) 
(14) 

-7/2 - 1 + 
(15) 



we obtain the desired result, i.e. that (|T2"]) coincides with the first of eqs. @. The numerical 
value is 

^(0) = 0.338096. (16) 

A second particular example is the following. For the case of a rectangle (of sides a and 
b) with Dirichlet boundary conditions, the spectrum of the Laplacian is X mn = 7T 2 (m 2 /a 2 + 
n 2 /b 2 ), and the zeta function 



Cr 



71 



-2s 



2^ I ~ + 52 



rrt,n=l 



H9 2 ' c(2s)+ ^& 



2s-l 



r(a-l/2) 
T(s) 



C(2a-1) 



(17) 



where, again, the corresponding asymptotic expansion for the Epstein zeta function in 
has been used. Taking now the derivative, we get 



C c (0) = ~ log(26) + ^ + £ K_ 1/2 (2nmna/b), 

z izo v m n=1 v m 



(18) 



which is best for numerical computations when a> b. In the particular case a = b (square), 
this reduces to 



77- / ji 

C(0) = -log(2a) + - + 2 £ J-K. 1/2 (2 

1 - m,n=l v 



mmn) 



which is just another expression for the same result obtained in |T(J (cf. eq. (All.3)) 

C,(0) = i log(2a) + i log(8vr) + ~ log 
The numerical value is, in both cases, 



(19) 



(20) 



CM = - log(2a) + 0.263672. 



(21) 



Notice, however, that for the general rectangle, expression flI8D is of much more practical 
use than the well-known one in terms of Dedekind's modular form 7], i.e. [ID, |TT] 

J_ 1/4 



C ec (0) = -log(a6)-log 



V2 V' 



r/(g) 



, v(q) = q 1/2i U(i-q m ), g = ex P I -2^-1 . 

(22) 

In fact, a few first terms of the series in ( |I8D suffice to obtain extremely accurate numerical 
results (just as for the case of the square). 



3 The vacuum energy density for Riemann surfaces in 
different dimensions 

We shall here calculate the vacuum energy density corresponding to a massless scalar field 
living in a Riemann sphere (manifold without boundary) or in a part thereof (namely a half- 
sphere, with Dirichlet or Neumann boundary conditions). The calculations become rather 
simple and precise numerical results are easy to obtain, by making exhaustive use of the 
formulas and considerations of the preceding section. But rather than attributing it to the 
specific manipulations we have carried out there (nothing special, in fact), this simplicity is 
to be interpreted as a success of the zeta-function regularization procedure itself (4|. 
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3.1 Dimension d=l 



This case is rather trivial and deserves no comment. For the one-dimensional Riemann 
sphere (no boundary, hence no boundary conditions), the vacuum energy density (or Casimir 
energy) is 

^ = 4 £»/S» = £ £>*«_!), (23) 

Z ' n=l ' n=l ' 

being r (here and in what follows) the radius of the Riemann sphere, and 

oo 

C(a) = E"-'. (24) 

n=l 

just the ordinary Riemann zeta function; therefore 

h , , 

Ex = -. 25 

12r 2 v ; 

For the semicircumference, the eigenmodes are shared by the cases of Dirichlet and 
Neumann boundary conditions, respectively. We obtain 

E° = E» = \. (26) 

i i 24r 2 V ; 

3.2 Dimension d=2 

For the ordinary Riemann sphere (i.e., the one which appears in string theory, see for instance 



10, 15]), we have, in the case of the half-sphere with Dirichlet and Neumann boundary 



conditions, 

£f = -AX>fA„ < = ~^X>f%, (27) 

i i 

respectively, where the the eigenvalues \ and eigenmode multiplicities mf, mf are 



\= V ^ r +1) , mf = Z, mf = / + l. (28) 
We get 

* oo oo 

E?> N = — 3 ( 2 D ' N (-l/2), ( 2 D (s) = + 1)]-, ( 2 N (s) = £(i + + 1)]-, (29) 

and by using our method as described in the preceding section, 

tfM = 2- + g fc , r J_^ fe) [C(2 a + fc -i)-H, 

6"W - *- + g arg^ K(2. + *-i)-i]. Po) 
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Now, when evaluating C2 ,N ( S = ~~ 1/2), a pole for k = 3 appears, in each case. The usual 
prescription in zeta function regularization is to take the principal part of the pole 0, |16, [17 . 
With this in mind, it is easy to obtain 

(— 1/2) = — — — — — + 0.033532, ^(-1/2) = — rr - 0.298630. (31) 

S2 v 1 ' 32(s+l/2) 2 v 1 ' 32(s + l/2) v ; 



This yields for the energy 



EP = 0.016766 - 4, E 2 = -0.149314 • 4- (32) 

iyO " /yd 



For the whole, ordinary Riemann sphere (no boundary), we have 

oo 

C 2 (s) = £(2/ + + 1)]"' = (ft*) + C*00 = -0.265096, (33) 
i=i 

and thus we obtain 

E 2 = -0.132548 • — . (34) 

Summing up, we see that in 2 + 1 spacetime dimensions, where space is the ordinary 
Riemann sphere (no boundaries), the Casimir energy density is negative. It is the sum of the 
energies corresponding to two half-spheres, one with Dirichlet and the other with Neumann 
boundary conditions on the one-dimensional boundary. The poles have opposite sign, and 
they anihilate when performing the sum. This is certainly an interesting result. It should be 
compared with the Casimir energy in 3 + 1 spacetime and boundary conditions imposed on 
a two-dimensional spherical surface. It is clear that the sign of the Casimir energy density 
for Riemannian manifolds can have important global cosmological consequences in plausible 
models of our universe. In fact, it could account for a lessening of the expansion ratio of the 
universe. 



3.3 Dimension d=3 

The three-dimensional Riemann sphere is a manifold without boundary that could perfectly 
well correspond to the spatial part of our universe, as a whole. The eigenvalues of the 
Laplacian operator are A 2 = 1(1 + 2)/r 4 , with degeneracies rrii = (I + l) 2 . Thus, the vacuum 
energy density for a massless scalar field is given by 

E 3 = ~ Us = -1/2), CsOO = EC + l ?W + 2)]- s . (35) 
/r i=i 
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and 



We can write, 

COO = E ^ s) d - r 2 r s = E Lrn I ^ 2g + 2k ~ 2 ) - ^ ^ 

;=2 fc=o K!i \ L ~ s ~ K > 

C3(-l/2) = g y^[ C (2* - 3) - 1] = - Ii ( JTT ^ - 0.411502, (37) 
which has a pole, for k = 2. Doing as above, we obtain the following numerical result: 

E 3 = -0.205751 • -j. (38) 

3.4 Dimension d=4 

For the four-dimensional Riemann sphere, the corresponding eigenvalues and multiplicities 
are Xf = 1(1 + 3)/r 6 and m ; = (I + 1)(Z + 2)(2/ + 3)/6. The vacuum energy density is now 

h 1 00 

E 4 = -— Us = -1/2), Us) = - E(^ + m + 2)(2/ + 3)[I(Z + 3)]"'. (39) 



We can write, 



1 00 1 00 f— ll^rfl — s 1 ) 

W .) = |g^_ 1/ 4)(u'-»/4r = igL|Sl-g 

9 x fc 



x (^-J [C(2s + 2A; - 3, 5/2) - C(2s + 2* - 1, 5/2)/4], (40) 
being u — I + 3/2 and £(s, a) Hurwitz's zeta function 

oo 

C(a,o) = 5> + a)-'. (41) 

n=0 

Again, nothing else has been done here but to apply the procedure as described in the 
preceding section. We thus obtain 

C4-1/2) = \ E (D k lC{2k " 4 ' V2) ~ ~ 2 ' V2)/4] = -°- 424550 - ( 42 ) 

It might seem that the term (9/4) fc could spoil convergence. This is not true: as always, we 
get a quickly convergent series. This is guaranteed in advance by the procedure itself, but 
it is rewarding to check this property numerically and see explicitly the rapid convergence 
(warning to the reader: one cannot use Mathematica to compute this expression since the 
Hurwitz zeta functions ((s,a) are very ill defined in this program for values of a > 1; this 



applies at least to version 2.0). Another important surprise is the fact that expression (4"2|) 
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is finite, no pole appears in this case (contrary to previous situations, in (|42|) there is no 
principal part reduction). 
Finally, 

E 4 = -0.212275 ■ 4- (43) 

Before closing this section, the following observation is in order. As explained before, an 
alternative treatment of the zeta functions above would be simply to split the polynomial 
in powers of the summation indices and then use the method of |§ . It is easy to check that 
this procedure is much more lengthy than the one developed here. On the other hand, the 
cancellation of poles in this method must be done explicitly (resorting to expansions around 
all poles and zeros), while it is immediate in our procedure (actually, no pole is ever formed). 
We conclude that the most direct way (the one we use) turns out to be here, at the same 
time, the shortest, most rigorous and best suited for numerical evaluation. 

The numerical results corresponding to the different cases here considered are depicted 
in Fig. 1. The vacuum energy density for a massless scalar field living in Riemann spheres 
(no-boundary manifolds) is represented as a function of the space dimension, d, in units of 
h r~( d+1 \ for d = 1,2,3,4. We just observe that the function is monotonically decreasing, 
showing no other remarkable feature than the small plateau for d = 3,4. In fact, the values 
for the three- and four-dimensional Riemann spheres are very similar. There is a sort of 
stabilization of the Casimir force for this number of dimensions. 

4 The vacuum energy density for cylinders of Riemann 
surface section 

As a sort of simplified study of stability of the vacuum energy density against deformations 
of the space manifold considered, we shall now investigate how the preceding values change 
when we consider the cylinders, with and without boundary conditions, which are obtained 
by adding a flat suplementary dimension to the Riemann spheres and half-spheres considered 
above. 

4.1 Dimension d=2 

We start with the case of a two-dimensional cylinder whose sections are one-dimensional Rie- 
mann spheres (no boundaries) or semicircumferences with Dirichlet and Neumann boundary 

11 



conditions, respectively. For the two last cases, we have (for a comprehensible treatment, 
see pi) 

E ™ - ^ C * £ - ^ = (44) 

being 



n=l,0 



CxT («) = £ (- 2 ) 1/2 ^ = £i ^ 1 C(2, - 1). (45) 

1 l s J n =l,Q 1 \ S ) 



We obtain immediately (no pole appears): 



E D,N = Rg(2) = _ Q 0Q24 ft = _ Q 004g ft 



4.2 Dimension d=3 

For a cylinder made up of ordinary Riemann half-spheres with Dirichlet or Neumann bound- 
ary conditions, what we have is 



being 



Cg 1 (^) = r(g rr ' /2) f:W + l)] 1/2 ^ dTiW = E(^ + l)[^ + l)] 1/2 - S - (48) 

We obtain 

<£(.) - ^{^ + f ^^^Kp.,*-,-!,}, 
c?l(s) = E^iZH) j 2 3/ 2 -. + 1 ^ff_-.) t) [c(2 . + 1 _ 2) _ uj . (49) 

Here a pole appears in both cases, with the same residue — 1/(60a/7t). Again, one must be 
here very careful when putting s = — 1/2 in these expressions, since one has to expand the 
zeroes and the poles of the gamma and zeta functions in terms of s + 1/2. After doing so, 
and using the following table of values for the derivative of the zeta function |20| 



we obtain 



C'(-l) = -0.16542115, C'(~2) = -0.0304, £'(-3) = 0.0054, 

('(-4) = 0.0080, C'(-5) = 0.00066, . . . (50) 

E° = -0.0042 • 4, £fi = -0.0073 -— A , E 21 = -0.0115 ■ 4- (51) 
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4.3 Dimension d=4 



For the case of the four-dimensional cylinder whose sections are three-dimensional Riemann 
spheres, we will only consider the no-boundary situation. Now 

£3,1 = ^ /_ dk g(Z + l)VA: 2 + /(/ + 2) = ^ CsAs = -1/2), (52) 

being 

- ^ £ KM- - ^ £ g^[C(— 3)-]- 

(53) 

This case is also finite, no pole appears when computing £ 3j i(— 1/2). On the other hand, it 
turns out that the final value for the energy density is remarkably small: 

£3,1 = -0.0115 • A (54) 

actually, it exactly coincides with the value obtained in the former case. There seem to be 
an intriguing stability of the vacuum energy density for these manifolds at this number of 
dimensions, more remarkable because the formulas for the zeta functions leading to these 
values look very different indeed (actually, the three-dimensional case needed a principal 
part evaluation of the pole while the four-dimensional case is finite) . Such a plateau almost 
appeared for Riemann spheres, but here the values obtained are identical and the plateau is 
completely horizontal (see Fig. 2). 

4.4 Dimension d=5 

Finally, we shall consider the five dimensional cylinder whose sections are four dimensional 
Riemann spheres. Here 

where 

C4,i(*) = r( *~A /2) f> 2(1 -V - V4)[l - 9/(4u 2 )]V2- u = I + 3/2. (56) 
61 \ s ) 1=1 

Using the same methods as before, we obtain, 

r(s-l/2) ~ (-l)T(3/2-s) f^ k 



(57) 
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Here we get again the usual pole. The final value is 

En = -0.0218- — . (58) 

The results of this section are summarized in Fig. 2, where the vacuum energy density 
corresponding to a massless scalar field living in cylinders whose sections are Riemann spheres 
(again manifolds without boundary) is represented as a function of the space dimension d 
and in units of % r~^ d+l \ for d = 2,3,4,5. A similar monotonic behavior with dimension 
as in the case of spheres is observed. Moreover, we notice an absolute stabilization of the 
numerical value of the energy density, which is exactly the same for d = 4 as for d — 3 
(sections are two- and three-dimensional Riemann spheres). Also, the case d — 4 is somehow 
special in the sense that it gives a finite result, no pole arises there. To be remarked as well 
is the fact that the Casimir energies for these cylinders are almost exactly by one order of 
magnitude smaller (in absolute value) than the ones obtained for the corresponding Riemann 
spheres. 

5 Discussion and conclusions 

We have calculated in this paper the vacuum energy density (also called Casimir energy 
density 0) corresponding to a massless scalar quantum field living in different no-boundary 
universes, as a function of the number of dimensions. For models of the universe we have 
chosen mainly manifolds without boundary, since they seem to be the most accepted nowa- 
days and, among them, Riemann spheres, since they are to be counted among the simplest 
and most important of these manifolds. However, several manifolds with a boundary, with 
Dirichlet and Neumann conditions, respectively, have also been considered. We have studied 
the variation of the vacuum energy density with the dimension of the space, hoping to find 
some characteristic that may singularize some of the manifolds considered or a particular 
dimension, among all of them. For this purpose, we have also compared the results obtained 
for spheres with the ones corresponding to cylinders of spherical section. It is in this respect 
that an alternative title of our work could have been: search for the most stable quantum 
universe without boundaries (stability in the sense of the vacuum energy density). 

As for the methods of analysis employed, they have consisted in zeta function regular- 
ization supplemented with a simple and natural binomial expansion, which has been shown 
to be both rigorous and very well suited for performing the numerical calculations. The 
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principal part prescription has been used to deal with the poles that appeared in several 
cases. 

The most remarkable results of our investigation are the following, (i) The absence of 
poles (completely finite result) for the case of the four- dimensional Riemann sphere and also 
for the four-dimensional cylinder of three-dimensional Riemann spherical section. Also in 
the low dimensional cases (d = 1, 2) poles are absent, but this is not true for the intermediate 
case d = 3. Of course also in this case the principal part prescription gives a well-defined, 
finite result, but this does not seem to be quite as satisfactory (for the very purists), (ii) The 
exact coincidence of the results corresponding to a 3d and a 4d cylinder even if the initial 
expressions are completely different, and the first result is obtained only after subtracting 
the pole, (iii) The fact that the vacuum energy density for cylinders is by an order of 
magnitude smaller than the one corresponding to spheres of the same dimension. Since all 
the results are negative, this means that when the no-boundary manifold is not spherical, the 
attractive Casimir force diminishes considerably, (iv) Finally, the elegance, simplicity and 
mathematical rigour of our method, as compared with other approaches, is very remarkable. 
As discussed before, the merit for this is to be given in full to the procedure of zeta function 
regularization itself, which is being confirmed as the most beautiful and useful existing 



regularization procedure [21 . 
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Figure captions 



Figure 1. The vacuum energy density for Riemann spheres (no-boundary manifolds) as a 
function of the space dimension, d, in units of h r~( d+1 ), for d = 1,2,3,4. The function is 
monotonically decreasing, showing no other remarkable feature than the small plateau for 
d — 3,4. In fact, the values for the three- and four-dimensional Riemann spheres are very 
similar. There is a sort of stabilization of the Casimir force for this number of dimensions. 

Figure 2. The vacuum energy density corresponding to cylinders whose sections are Rie- 
mann spheres (again manifolds without boundary), as a function of the space dimension d 
and in units of h r~( d+1 \ for d = 1,2, 3, 4. We see the same monotonical behavior as in the 
case of the spheres but it is remarkable the clear stabilization of the numerical value of the 
energy density which is exactly the same for d — 4 as for d = 3 (cylinders whose sections are 
two- and three-dimensional Riemann spheres, respectively). 
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